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The influence of Raman scattering and higher order dispersions on solitons and 
frequency comb generation in silica microring resonators is investigated. The 
Raman effect introduces a threshold value in the resonator quality factor above 
which the frequency locked solitons can not exist and, instead, a rich dynamics 
characterized by generation of self-frequency shifting solitons and dispersive 
waves is observed. A mechanism of broadening of the Cherenkov radiation 
through Hopf instability of the frequency locked solitons is also reported. 
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1. Introduction 

Generation of broad frequency combs in microring res¬ 
onators provides an attractive tool for many practical 
applications, where a miniature source of broadband reg¬ 
ularly spaced spectral lines is required [1]. Typical exper¬ 
iments on microresonator combs have been performed 
with silicon nitride and silica glass rings and spheres 
EHU- Four wave mixing and soliton formation are the 
most important building blocks of microresonator combs, 
which makes an obvious link to the supercontinuum gen¬ 
eration widely studied in photonic crystal fibers 0 0 . 
Comb generation has been optimized, e.g., by engineering 
the group velocity dispersion (GVD) in order to tune the 
modulational instability (MI) gain 0 7, 8] and by pump¬ 
ing close to the zero GVD points to benefit from efficient 
Cherenkov radiation by solitons n mu]. Note, that 
a key ingredient for supercontinuum generation in silica 
glass fibers is the interplay between the self-frequency 
shifting Raman solitons (SFSR solitons) and dispersive 
waves 0. Raman effect in the microresonator comb gen¬ 
eration has received relatively small theoretical attention 
so far mm- Raman scattering in silica-based micro¬ 
cavities has been exploited to obtain Raman lasing, see, 
e.g., mm- However, detailed understanding of its role 
appears to be important from the fundamental point of 
view and also to promote use of promising silica based mi¬ 
croresonators in ongoing research dM3- We also note 
that mode-locking in active laser oscillators with Raman 
nonlinearity has been recently studied using generalized 
nonlinear Ginzburg-Landau equations Hoi ED- 

In this work we focus on the comb generation through 
the soliton formation and analyze the existence and sta¬ 
bility of cavity solitons in the presence of the broad band 


Raman gain and higher order dispersion (HOD) effects. 
Whereas in optical fibers a cw pump in the anomalous 
GVD range (and close to the zero GVD point) automati¬ 
cally yields supercontinuum generation driven by soliton 
dynamics [3 E>], the role of solitons in the microresonator 
combs is more subtle, since the soliton existence and sta¬ 
bility in resonators is not universal and critically depends 
on parameters. Trains of solitons associated with micro¬ 
resonator combs have been previously studied without 
the Raman effect in E2H2I] • Though, even in this case 
the impact of soliton bifurcations on the comb regimes 
still has not been explored fully. Below we demonstrate 
that the comb formation in silica micro-resonators with 
quality factors Q order of 10 6 is strongly affected by a 
type of Raman solitons, which frequency is locked to a 
specific value due to strict balance between the Raman 
gain and cavity loss. We term these solitons as frequency 
locked Raman (FLR) solitons. While in the silica micro¬ 
resonators with Q rsj 10 8 and above [T7HT9] . the low cav¬ 
ity loss is not able to balance Raman gain and no exact 
soliton solution has been found, while the comb forma¬ 
tion is dominated by the self-frequency shifting Raman 
(SFSR) quasi-solitons, similar to the ones observed in free 
propagation experiments in optical fibers 0 . Both types 
of solitons can emit resonant radiation contributing into 
shaping of frequency combs. 


2. Lugiato-Lefever equation and Raman effect in 
microring resonators 


* Corresponding author: carles.milian@cpht.polytechnique.fr 


To describe light dynamics in silica microrings, we use 
the equation for the envelope E of the intracavity field 
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Ee l P° z + c.c., where /3q = /3(c<Jo), 
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Here = dpiu = (d n uj/df3 n ) OJ=OJoJ (3 is the propagation 
constant of the waveguide the resonator is made from, uo o 
is the reference cavity resonance and the other resonances 
are given by u q = ujq+^u' q 2 +y^uj'"q^^r ..., 
where q is the modal number q = 0, ±1, ±2, ±3,... 
counted relative to the pump at q = 0 . a/ = c/n g , 
c j" = — (c 3 /n^)9 w (n^/c), where n g is the group index 
at co>o and c is the vacuum speed of light. T is the rate 
of photon loss from the cavity and A is the dimension¬ 
less amplitude of the pump field, r is the pump coupling 
rate and u p = 2ttc/\ p is the pump frequency, which is 
detuned from ujq roughly within the half of the cavity 
free spectral range (FSR). \E \ 2 is the dimensionless in¬ 
tracavity intensity and g is the dimensionless nonlinear 
parameter, see [251, 26] for scaling, t is time and z is 
the coordinate along the cavity varying from 0 to 2ttR, 
where R is the radius: E(z = 0) = E(z = 2ttR). x is the 
Raman part of the nonlinear susceptibility with the stan¬ 
dard parameters used for silica glass: /i = 0.18, 7 # = (16 
fs) -1 , Q r = [y R /A + W fl] 1/2 J and W -R = ( 12 - 2 fs) _1 - 
Equation (1) has been derived in details in [26] in the 
context of ring microresonators. In 2_6 it was pointed 
out that this equation can be straightforwardly gener¬ 
alised to include other physical effects and, in particu¬ 
lar, Raman effect. We have followed the approach in 
Ref. [26] with the only difference that the nonlinear po¬ 
larisation includes in our case the Raman term obeying 
the standard Raman oscillator equation E3- Indeed, the 
Raman part of the material susceptibility in the classi¬ 
cal approximation is well known to be described by the 
oscillator equation driven by the intensity of the applied 
field [271 EE], which is our Eq. (2). A solution to the 
Eq. (2) can be expressed in an integral form and it has 
been a preferred method of describing Raman effect in su¬ 
percontinuum generation in optical fibers 0, cf. [28]. A 
model similar to ours in a multi-mode coherently pumped 
Raman laser has been previously reported in [29] , where, 
however, the Raman effect was treated quantum mechan¬ 
ically using a 3-level model and dispersion of the cavity 
modes was assumed linear. 

Note here, that there are two approaches widely used 
in literature to describe an optical field evolution in res¬ 
onators. One is using time as an evolution variable 
and relies on the approximation of the modal disper¬ 
sion through the expansion of the frequency into series in 
the modal number (propagation constant) and the other 
one is using the distance along the cavity length as an 
evolution variable and hence relies on the expansion of 
the modal number as a function of frequency. Recently 


the first approaches was used by Chembo and co-author 
[25l 126] , by Kippenberg and co-authors nang and by us 
m, while the second approach is most often associated 
with a series of papers by Haelterman and co-authors 
and by many others In the lat¬ 

ter case, periodic boundary conditions still have to be 
applied in the spatial coordinate, which, however, enters 
in the model through the first derivative. A significant 
advantage of the approach adopted here is that it allows 
a straightforward theoretical study of instabilities, since 
the instability develops in time and periodic boundary 
conditions in space are easily implemented. While an at¬ 
tempt to study how instability develops along the cavity 
length, which is a periodic coordinate, though not im¬ 
possible, runs into unnecessary mathematical difficulties 
1321 . Also note, that we avoid using here the terminology 
of the ’slow time’ (which is the distance along the cavity 
scaled to the group velocity) and ’fast time’ (which is the 
physical time) [31, . 

Introducing normalization E = -^==^/e~ luJpt ^ x = 
t = Tt with r = [27r R]/[c/n g \ (so that the 
group velocity coefficient becomes unity), and distance 
z = ZL with L = 2i tR (R is the resonator radius) we 
get a more handy form of the dimensionless generalized 
Lugiato-Lefever [33] equation (gLLE) 

M 

id T v - £ B m (-id z ) m v + (3) 

m= 1 

+(»7 - <5 + 2[(1 - / x)|^| 2 + W ])¥ + h = 0 

^d 2 T w + ^d T w + {w-m 2 ) = ^ (4) 

R R 

where B rn = /{27iR) rrl /m\ ( B\ = 1), 7 = Tr. S = 

(cjo — cj p )t, h = VTyJgujQrA. Note, that periodic bound¬ 
ary conditions along the resonator circumference make 
it natural to use the model with time, not space, being 
an evolution coordinate and to treat Raman polarization 
through the separate equation Eq. (4) and not the con¬ 
volution integral as in fiber propagation problems. For a 
silica microresonator with R = 400 gm and the waveg¬ 
uide diameter 1.2 fi m considered below and pumped at 
A p = 1.47 \i m, we have: FSR is c/(n g R) = 480 GHz 
or in dimensionless units Sfsr = 27r, round trip time 
r = 13.12 ps, and quality factor Q = t(2ttc/ Ap/ 7 ), which 
we vary from 10 6 upto 10 8 . 

3. Frequency-locked and self-frequency shifting 

Raman solitons 


Exact soliton solutions of Eqs. (3)-(4) are sought in the 
form 4/(T, Z) = ^(x), W(T,Z) = w(x) where x = Z — 
(v + 1 )T and w obey: 


- ivd x ip + B 2 dlil} - iB 3 dlip + 

(*7 — S + 2[(1 - n)\4>\ 2 + w\)ip + h = 0, 

{y + 1) 2 2 (1 + u ) 7 h | | 2 

002 d x w - ~Fa - d x w + w = n\i>\ , 


T 2 n 2 R 


tQ, 2 r 


(5) 

(6) 
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Fig. 1. (Color online), (a) soliton branches (B 3 = 0), with and without Raman, and bistability for S = 0.1, 7 = 0.02 
(Q = 8.5 x 10 5 ), and B 2 = 2.38 x 10 -7 . Thick traces denote soliton stability, (b) spatial and (c) spectral representation 
of the solitons with h = 0.0045 (gray) and h = 0.0089 (black). Inset in (c) shows several cavity modes around the 
pump wavelength, (d) Hopf instability growth rate for the solitons in (a). Inset shows Raman soliton and Hopf 
modes spectra (linear scale) for h = 0.0089. (e,f) Spatial soliton propagation for h = 0.0045 (stable) and h = 0.0089 
(unstable), respectively. Dashed lines mark the soliton velocities calculated with the Newton method, demonstrating 
exact coincidence. 



Fig. 2 . (Color online). Soliton momenta branches vs loss 
with /i = 0.18, £>3 = £>4 = ... = 0 for single (black) 
and double (orange, light) peak solutions. Solid thick 
lines show stability regions and ”H” marks the onset for 
Hopf instability. Dashed vertical line at low 7 marks 
the threshold quality factor separating the unstable exis¬ 
tence (see Fig. 3(b)) and non-existence, Figs. 3(c)-3(d), 
regions. ”T” indicate turning points (local minima of 7 ) 
that are not further explored in this work. Insets show 
real (dashed) and imaginary (solid) spatial profiles of soli- 
tons around the double turning point: A—)► Bl, B 2 . 


where, for the sake of simplicity, we are restricting our¬ 
selves to the third order dispersion (TOD). In what fol¬ 
lows we rely on numerical solutions of Eqs. (5,6) for the 
soliton profiles found by a modified Newton method and 
on the numerical stability analysis based on finding eigen¬ 
values and eigenvectors of the Jacobian derived from the 
linearization of Eqs. (5,6) around the soliton solution. 

The amplitude of the single mode state (flat solution) 
at the pump frequency, found from (iy — S + 2|'0o| 2 )V ; o + 
h = 0 , is multivalued (bistable) in the soliton existence 
range and the solitons are nested on the background given 
by the root with the smallest value of IV’ol, see Figs. 
1(a) and 1(b). Figure 1(a) shows bistability and soli- 
ton branches found with (i = 0 (upper) and fi = 0.18 
(lower) and for 7 = 0 . 02 , S = 0 . 1 , B% = 0 . These soli¬ 
tons exist for some special values of v fixed by the choice 
of the other parameters and found self-consistently with 
the soliton profile, i.e., they are the co-dimension 1 struc¬ 
tures in the bifurcation theory terminology. Physically, 
different v’s imply different carrier frequency of the soli¬ 
tons and hence we term such solitons as the frequency- 
locked Raman (FLR) solitons. A transition between sta¬ 
ble and unstable regimes for FLR solitons is shown in 
Fig. 1(d). The instability is the oscillatory (Hopf) one, 
as is well known in the cavity soliton context, see [34] 
and references therein. Figures 1(b) and 1(c) show spa- 
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Fig. 3. (Color online). Intracavity field after seeding with a pulse plus: fi = 0.18, h = 0.003, 8 = 0.1, B 2 = 2.38 x 10 -7 , 
HOD= 0 . (a-c) space-time evolution for 7 = 0 . 012 , 7 = 0.005, 7 = 0.0001 (or Q = 1.4 x 10 6 , Q = 3.3 x 10 6 , 
Q = 1.7 x 10 8 ), corresponding to the excitation of the stable and Hopf unstable FLR solitons, and to the SFSR 
solitons, respectively, (d) spectral representation of (c). (i-iv) Spatial profiles corresponding to the dashed lines in 
(c), illustrating the growth of the NLS-like accelerating solitons. The field modulus, |^|, (black) and its imaginary 
part, Im(^), (red, lighter) are plotted. Positive (negative) values of Im(^) are marked by the dark blue (light gray) 
areas as to illustrate where localized amplification (absorption) is possible. 


tial and spectral profiles of FLR solitons. Raman effect 
predictably induces the red-shift of the carrier frequency 
corresponding to v < 0, cf. Fig. 1(c). Propagation dy¬ 
namics in space-time is shown in Figs. 1(e) and 1(f) for 
solitons with h=0.0045 (stable) and 0.0089 (Hopf unsta¬ 
ble). Persistent oscillations of solitons are also univer¬ 
sally observed closer to the Hopf threshold [30, 35, 36] . 
see below. 

FLR solitons can be found only when losses are present 
(7 > 0) and exceed some critical level, see Fig. 2, while 
for fi = 0 solitons in the Lugiato-Lefever model exist 
all the way down to 7 = 0 HUE]. This can be seen 
from the momentum, M m\ Jq {'ifrdx'ip* — c.cjdx, balance 
equation, which for the time-independent traveling wave 
forms reads as 

7 M s = - [ wd x \ip s \ 2 dx. (7) 

Jo 

Thus localized solutions with w 7 ^ 0 give a finite right- 
hand side and hence can only exist above some thresh¬ 
old in 7 . Figure 2 shows single and double peak soliton 
branches in the M, 7 plane in the presence of Raman ef¬ 
fect and without HODs, for h = 0.003, 5 = 0.1. While 
the single hump branches together with the double hump 
ones A and B 2 are forming a well known snaking bifurca¬ 
tion diagram, the double-hump branch HI is likely to be 
a part of a more complicated ’snakes-and-ladders’ struc¬ 
ture, see, e.g., [38] . 

Figures 3(a)- 3(d) show propagation dynamics seeded 


by a pulse for three different values of loss correspond¬ 
ing to the stable FLR soliton, Fig. 3(a), Hopf unstable 
FLR soliton, Fig. 3(b), and to the low loss (high Q) 
range, where FLR solitons do not exist, Figs. 3(c) and 
3(d), respectively. The first two cases demonstrate con¬ 
vergence of the input pulse to either a stable soliton, 
Fig. 3(a), or a Hopf unstable one, Fig. 3(b), leaving 
the cavity in the single mode state. Exciting the cavity 
with Q > Q t h = 5.2 x 10 7 , Figs. 3(c) and 3(d), ex¬ 
hibits a very different dynamics. In this case we have 
observed cascaded generation of localized pulses experi¬ 
encing continuous frequency shifts away from the pump 
field towards longer wavelengths. Thus, in the high Q 
regime, instead of the FLR solitons with v = const , 
we have self-frequency shifting Raman (SFSR) quasi- 
solitons, which are essentially solitons of the nonlinear 
Schrodinger (NLS) equations with complex dynamics due 
to interaction with pump and other dispersive waves. 
Qualitatively, it is not surprising that FLR solitons cease 
to exist in the limit of low loss, pump, and detuning, i.e., 
as Eqs. ( 1 , 2 ) become the NLS equation, which has SFSR 
soliton solutions [ 6 ]. It is, however, not only when the 
propagation is close to that in fibers that the accelerating 
solitons appear. As discussed above (and in agreement 
with our results, see, e.g., Figs. 3(c) and 3(d)), for 7 —>> 0 
no FLR solitons can exist, regardless the strength of the 
pump or the magnitude of the detuning. 

While the dissipative solitons have been observed in 
silica fiber loops Isa am, the influence of Raman effect 
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Fig. 4. (Color online). FLR soliton profiles in spatial (a-c) and (d-f) spectral domains for /i = 0.18, 7 = 0.02, 5 = 0.1, 
h = 0.007. (a,b) is for the positive TOD, B 3 = 0.8 x 10 -10 , giving the quiescent soliton, and (c,d) is for the negative 
TOD B 3 = -2.48 x 10 10 . Insets in (a,b) show the Raman oscillator field, w. Solid lines in (d,e) correspond to spectra 
without Raman effect. RRq 2 ) mark the resonant radiation (RR) roots. FWM peaks resulting from the interaction of 
pump and RRi are also labeled. (c,f), B 3 = 2.48 x 10 —10 , correspond to the oscillatory soliton in Fig. 6 at T = 920. 
Vertical solid lines mark the predicted resonances and the dashed the zero GVD. Numbers in (f) correspond to m in 
Eq. (9). 


on them has not been reported so far, which can be due 
to pico-second pulse durations and low quality factors 
resulting in frequency-locked propagation with negligible 
frequency offsets. Note here, that the frequency locked 
soliton states in the SiN ring resonator have recently been 
modeled numerically m using an approach when the 
first derivative in gLLE is taken in the propagation coor¬ 
dinate. The Raman effect was also included there, how¬ 
ever, its impact on the soliton dynamics was found to 
be only quantitative and did not receive an attention be¬ 
yond a brief comment. This is likely due to the fact that 
the Raman line in SiN is very narrow if compared to the 
silica one, so that its impact on the soliton spectrum and 
dynamics is minimal. 


Note, that the cascaded SFSR soliton generation in 
Fig. 3(c) is accompanied by the increase of the power 


integral, P = \ij;\ 2 dx, which evolution is given by: 

\d T P = - 7 P + h f Im(ip)dx. ( 8 ) 

2 J 0 

Figures 3(i)-3(iv) show spatial field profiles at distances 
marked by (i)-(iv) in Fig. 3(c). Generation of pulses 
that are seen to detach from the input one, see Figs. 3(i)- 
3(ii), is accompanied by the growth of the power integral, 
which is possible only if, at least locally, f Im{^)dx > 
r yP/h . The latter condition is well satisfied for the grow¬ 
ing pulses, which imaginary parts are shown in Figs. 
3(iii) and 3(iv). These keep the ^ ir/2 phase differ¬ 
ence with the pump (exhibiting a pronounced imaginary 
part) sustaining the growth condition. Note here that fre¬ 
quency locked combs were demonstrated in high-Q fused 
silica micro-rings in, e.g. 14M431 , their relation to solitons 
have not been investigated in those papers and additional 
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experimental and theoretical studies are required. 

4. Resonant radiation of stable and Hopf unsta¬ 
ble FLR solitons 

TOD gives rise to cavity solitons with radiation tails, 
which have been recently reported numerically mmm 
05] and observed experimentally da is]. Resonant ra¬ 
diation is expected in general for any localized nonlinear 
wave, including bright and dark solitons and shock waves 
m im ns]. Emission of the resonant radiation is known 
to induce the spectral recoil on the soliton core shifting 
the soliton carrier frequency away from the zero GVD 
point [61 09]. Figure 4 shows typical spatial and spec¬ 
tral soliton profiles with /i, B% ^ 0. The two resonant 
radiation (RR) peaks are in agreement with the expres¬ 
sion for the resonance conditions derived in m- Other 
spectral features occurring because of the nonlinear four 
wave mixing (FWM), e.g., between the radiation and the 
pump, are labeled as ”FWM”. A comparison of the soli¬ 
ton spectra with (red area) and without (blue solid curve) 
Raman does not reveal significant changes in the solitons, 
apart from a slight tuning of the radiation frequency and 
amplitude due to the red-shift of the core. Raman effect 
enhances the radiation tail for L > 3 > 0 (/% < 0 ) since 
in this case the spectral maximum of the FLR soliton is 
shifted closer to the zero GVD wavelength, see Fig. 4(a) 
and 4(d). The opposite happens for L > 3 < 0 (/% > 0), 
similarly to what is known in fiber propagation [6 . An 
interesting feature of the FLR solitons is the existence 
of a quiescent soliton, v = 0, for B% > 0, see Fig. 5(a). 
This manifests an exact balance between the radiation 
induced spectral recoil and the Raman effect. 

Figure 5(b) compares the growth rate of the Hopf in¬ 
stability in presence and in the absence of the Raman 
effect as a function of L > 3 and for FLR solitons that are 
either stable (h = 0.0055) or Hopf unstable (h = 0.007) 
when L >3 = 0. For /i = 0.18, i.e., due to Raman effect, 
growth rate is not symmetric with respect to L> 3 . We 
also note that large values of TOD have unambiguous 
stabilizing influence on the solitons. Sufficiently close to 
the Hopf instability threshold the cavity solitons with 
and without Raman effect form stable breathing struc¬ 
tures oscillating with the Hopf frequency D#, see Figs. 
4(c), 4(f), and[ 6 j A notable feature here is that the res¬ 
onant radiation of the breather is a set of peaks form¬ 
ing the Cherenkov comb or the resonant radiation comb, 
see Figs. | 6 jc) andj 6 ^d). The Cherenkov comb is much 
broader than the single Cherenkov peak, thus transition 
from stable to Hopf unstable solitons significantly en¬ 
hances the wing of the overall resonator comb shaped by 
the dispersive wave effects, see Fig. HD and cf. Figs 4(d) 
and 4(f). These transformation of the Cherenkov spec¬ 
trum from a single line into a comb, cf. Figs. 4(a) and 
4(d), Figs. 4(c) and 4(f), can also be seen in Fig. 8 (b) 
in m- Time domain spectrum of the oscillating FLR 
soliton, see Fig. [ 6 jb), shows a number of well defined 
spectral components separated by Q#. Each of these 


components serves now as a source of the Cherenkov ra¬ 
diation, so that the radiation phase matching condition 
derived in m (see Eq. ( 6 ) there) for the carrier fre¬ 
quency of a stable cavity soliton has to be generalised to 
the multi-frequency structure, corresponding to the soli¬ 
ton breather. Dimensionless form of this generalization 
is readily expressed as: 

0 = mUfl - vQ m + BzQ^ ± (9) 

±\/i<5 + -B 2 <3m — 4|'E ' 0 | 2 ] 2 — 4|\l/o| 4 , m = 0, ±1, ±2,... 

Here mfln are the frequencies of the soliton breather 
components and Q m is the modal number correspond¬ 
ing to the resonance radiation, the associated resonance 
frequencies in physical units can be computed as ujq + 
Q m [v + l}c/[2irRn g \ [lT . An identical generalisation of 
the resonant condition for the multi-soliton breathers in 
optical fibers (nonlinear Schrodinger equation) has been 
introduced in [50;. Fig. 5(c) shows graphical solution 
of Eq. [lO] and Figs. 4(f), [ 6 ][d) are mapping the pre¬ 
dicted resonances onto the numerically observed ones. 
We note, that new resonances appear not only around 
the Cherenkov frequency, but also across the spectrum 
of the soliton core, cf. Figs. 5(c) and 4(f). Considering 
complexity of the dynamics impacting the instantaneous 
shape of the soliton spectrum, the matching of the reso¬ 
nances found from Eq. (9) with the numerical data can 
be transformed from qualitative to quantitative if a fit¬ 
ting parameter controlling shift of the central frequency 
is introduced into Eq. (9). Related, but notably dif¬ 
ferent in its realization and spectral manifestation, phe¬ 
nomenon of multiple dispersive wave resonances in fibers 
with dispersion modulation have been recently explored 
in [5DE2]. 

5. CW excitation of frequency combs 

Observation of the microring solitons has been originally 
reported by simply tuning the pump frequency across 
the bistability range of the resonator [22 . Following this 
approach, we show in Fig. [7][a) a bistability diagram for 
the single mode state and the families of the FLR solitons 
traced as a function of the cavity detuning S. Soliton sta¬ 
bility intervals shown by the thick lines shrink when L > 3 is 
changing from positive to negative, due to interplay with 
the Raman effect, consistent with Fig. 5(b). Though the 
Raman effect does not impact the range of existence in 
pump, it significantly narrows it in S, cf. Fig. 1 (a) and 
0 a). Excitation of the resonator with a cw pump re¬ 
sults in a variety of regimes ranging from spatio-temporal 
chaos to stable soliton propagation observed for various S 
|9l [ 221 . In order to investigate how these regimes can be 
altered by the Raman effect we have computed disper¬ 
sion profiles, see Fig. 7][b), for two micro-ring resonators 
made of silica strands with two different diameters of 
4 and 1.2 /im with the ring radius R = 400 fi m. The 
4 /im diameter resonator is pumped at 1 /im close to 
its only zero GVD wavelength, when TOD is negative, 
so that towards the longer wavelength from the pump 
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Fig. 5. (Color online). Soliton velocity (a) and the Hopf instability growth (b) with (fi = 0.18, light blue) and without 
(/i = 0, black) the Raman effect as a function of the third order dispersion parameter. The curves in (a,b) emerge 
from those in Fig. 1(d) for h = 0.0055 (dashed) and h = 0.007 (solid): S = 0.1, 7 = 0 . 02 . (c) Phase matching 
diagrams from Eq. (9) for m G [—5,5] (see Fig. 6 for parameters). 



Fig. 6 . (Color online), (a) Oscillations in time T of the spectral maximum of the Hopf unstable FLR soliton. 
(b) Fourier spectrum of the oscillations shown in (a), revealing multiple frequency components separated by the 
Hopf frequency H#. (c) Periodic temporal dynamics of the spectrum of the Hopf unstable FLR soliton core and of its 
radiation tails. Dashed vertical line shows zero GVD wavelength, (d) Spectrum of the radiation tail showing formation 
of a frequency comb. Vertical lines show prediction of the new resonances from simple theoretical considerations, see. 
Eq. (9). B 3 = 5 x IQ” 11 , S = 0.1, 7 = 0 . 02 , h = 0.007, [i = 0.18. 


the GVD is always anomalous. The 1.2 fim resonator is 
pumped close to 1.5 fi m, where TOD is positive and the 
zero GVD point and normal GVD range are towards the 
longer wavelengths. 


For relatively low quality factors (Q ~ 10 6 ) dynam¬ 
ics in both cavities is similar. For detunings outside the 
soliton existence range in the interval marked as ’chaos’ 
in Fig. 0a) and within the FLR soliton existence inter- 







































































6. Conclusions 



Fig. 7. (Color online), (a) Bistability vs detuning and 
soliton branches for B 2 = 2.38 x 10 -7 and 7 = 0.01 (Q = 
1.7 x 10 6 ): (black) B 3 = fi = 0, (dark red) B 3 = 0,/i = 
0.18, and (light grey) L > 3 = ±2.48 x 10 -10 , /jl = 0.18. 
(b) GVD for the two silica strands with cross section 
diameters, d = 1.2 fim and d = 4 /am (see Figs. [ 8 | 
[To| ). Numerical simulations presented below use the full 
dispersion profiles in (b). 


val, but close to its left boundary, the development of MI 
leads to chaotic spatio-temporal patterns shown in Figs. 
Up), life), go), and ||g). As the detuning is increased 
and FLR solitons become more stable, the MI of the cw 
pump leads to formation of the sparsely separated soli¬ 
tons, which form a nice frequency comb. FLR solitons 
emit dispersive waves, which broaden and reshape the 
combs. All radiation peaks in the exact soliton spectra 
seen in Fig. 4 can be traced in Fig. [ 8 ] as well and are 
even more pronounced due to the large number of radi¬ 
ating solitons. Note, that the broad MI spectrum in Fig. 
[ 8 |c) masks the strong Cherenkov peak (RRi) at 1600 nm. 


Figures [9] and [lO] show comb generation in the same 
cavity geometries and for the same pump wavelengths, 
but when the FLR solitons do not exist since we have 
increased quality factor by two orders of magnitude to 
Q ~ 10 8 . Now cw MI generates SFSR quasi-solitons, 
which are able to reach large frequency detunings in the 
anomalous GVD range in the 4 /im resonator without 
been appreciably attenuated and at the same time to gen¬ 
erate significant amount of dispersive radiation into the 
normal GVD range, see Figs. [9ja) and §b). Increasing 
the pump levels leads to formation of shorter and more 
intense solitons and hence to stronger SFSR solitons and 
broader combs. The frequency shift of the solitons in 
the 1.2 /1 m resonator is, however, quickly arrested by the 
normal GVD range extending towards the longer wave¬ 
lengths and antagonistic with the soliton existence, Fig. 
[To] In fact, resonant radiation in this case amplifies to 
very large amplitudes, so that after some time, T > 300, 


it wipes out the solitons from the resonator, see Fig. 10 


Thus, when FLR solitons are replaced by SFSR solitons 
due to increase in the resonator Q the sign of the third 
order dispersion is playing a crucial role in shaping of the 
frequency combs. 


In conclusion we have studied the role of the Raman effect 
in frequency comb generation in silica microresonators. 
We identified two different regimes of comb formation 
associated with two different types of solitons. One type 
is the frequency-locked Raman solitons existing in res¬ 
onators with Q ~ 10 6 and the other is the self-frequency 
shifting Raman quasi-solitons dominating comb forma¬ 
tion when Q ^ 10 8 . In the former case the comb is 
formed by the static soliton spectrum and the resonant 
radiation tails, while in the latter situation the comb is 
dynamic and is expanding in time, similarly to the expan¬ 
sion of supercontinuum spectra generated in optical fibers 
[ 6 ]. Factors limiting this expansion in microresonators 
require further investigation. We are also reporting sig¬ 
nificant broadening of the Cherenkov radiation spectrum 
through the formation of the resonant radiation comb, 
when FLR solitons become Hopf unstable. Further inves¬ 
tigation of various aspects of the effects reported above 
is warranted. 

Finally, we acknowledge very recent investigations 
about microring Raman solitons which are of particu¬ 
lar interest for the current work. The soliton red-shift 
has been measured experimentally in amorphous SiN in 
Ref. |53|, which has a broad band Raman gain, similar 
to silica, albeit with different natural frequency, £ 2 #, and 
damping, j R . 
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